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The many-body reciprocal theorem and swimmer hydrodynamics 
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We present a reinterpretation and extension of the reciprocal theorem for swimmers, extending 
its application from the motion of a single swimmer in an unbounded domain to the general setting, 
giving results for both swimmer interactions and general hydrodynamics. We illustrate the method 
for a squirmer near a planar surface, recovering standard literature results and extending them to 
a general squirming set, to motion in the presence of a ciliated surface, and expressions for the 
flow field throughout the domain. Finally, we present exact results for the hydrodynamics in two 
dimensions which shed light on the near-field behaviour. 


I. INTRODUCTION 


Microorganisms and self-propelled particles move and 
interact in a fluid environment dominated by the effects 
of viscosity w- The long-range nature of Stokes flows 
suggest that hydrodynamics will have a significant effect 
on microorganism motion and behaviour. Examples of 
this are seen in the circular motion of E. coli near bound¬ 
aries 00, attraction to surfaces 00, synchronisation 
of cilia and flagella 00, interactions and scattering of 
pairs of swimmers jlOlllll . and collective dynamics of sus¬ 
pensions |l2tll4] . The difficulty in determining swimmer 
hydrodynamics theoretically comes because the motion 
of the organism, which we wish to determine, constitutes 
a boundary condition. An elegant method for overcom¬ 
ing this was introduced by Stone & Samuel [TFj, based 
on the Lorentz reciprocal theorem, a general relationship 
between two Stokes flows in the same domain but with 
different boundary conditions. This approach was orig¬ 
inally developed as a way of determining the motion of 
a single swimmer, with a significant feature that it by¬ 
passes the need to calculate the full hydrodynamics, and 
in this capacity has become a standard tool in the active 
matter literature [TM2U1 . 

In fact, the reciprocal theorem may be used to deter¬ 
mine essentially all aspects of the hydrodynamics of any 
swimmer problem, a point apparently not emphasised in 
the literature. Here we describe how this may be done 
and provide illustrations for a wide range of problems in 
two and three dimensions. We recover standard results in 
three dimensions for the interaction of a swimmer with a 
planar wall m and extend to the case of a swimmer in¬ 
teracting with an active carpet of cilia, showing that the 
surface activity may cause trapping or enhance deflection 
depending on the direction of the surface flow. 

The full utility of the reciprocal theorem is that prob¬ 
lems in swimmer hydrodynamics may be solved easily if 
a conjugate Stokes drag solution is known. We illustrate 
this in two dimensions giving exact results for problems 
covering the interactions of two circular squirmers, or the 
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motion of a single particle either inside or outside a cir¬ 
cular domain, or swimming close to a planar boundary. 


II. THE RECIPROCAL THEOREM 

The reciprocal theorem [22] relates two solutions of 
the Stokes equations in the same domain, D , but with 
different boundary conditions. The solutions, (it, <x) and 
(■u,<x), are related by the integral relation 




n, 


(1) 


where fi is the unit normal pointing into the fluid, u is 
the velocity and <x is the stress. For a collection of swim¬ 
mers moving with no force or torque at translational and 
rotational speeds Ui and due to surface slip velocities 
u S i , this splits as a sum over the boundaries, Si, of each 
of the swimmers to give 


X [ui ■ + a • /: 



( 2 ) 


where <l?i and T* are the force and torque on the particle 
in a conjugate solution (u, a). By an appropriate choice 
of forces and torques in the conjugate problem each com¬ 
ponent of the swimmer motion can be obtained. The 
right-hand-side expresses this swimmer motion in terms 
of the boundary data, the slip velocities u si , with the 
stress tensor of a conjugate Stokes drag problem playing 
the role of an integration kernel. 

Although only a mild generalisation of the result given 
in [15] , eq. ([2]) is tremendously general and allows more or 
less all aspects of swimmer hydrodynamics to be solved 
for. Loosely, this is as expected. The flow is given by the 
boundary data; the role of the reciprocal theorem is to 
identify the appropriate integration kernel. 

We illustrate this first for the case of a swimmer close 
to a planar boundary. The integration kernel is the stress 
tensor ef for the Stokes drag of a particle in the half-space, 
for which an approximate expression can be given in 
terms of fundamental flow singularities and the method 
of images. Retaining only the leading order contribution 
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FIG. 1: The effect of a wall on swimmer trajectories. Close 
to a no-slip wall, a squirmer oriented parallel to the wall is 
repelled, following the black trajectory. This repulsion may 
be suppressed by a carpet of cilia on the wall generating a 
metachronal wave moving in the same direction, which can 
trap the swimmer in an advected orbit, shown by the red 
trajectory. If the metachronal wave opposes the swimming, 
the repulsion is enhanced, shown in blue. The black arrows 
represent the activity at t = 0, resulting in the flow field, 
which decays exponentially with distance, shown in grey. 


of a point force and torque in the half-space [23], the in¬ 
fluence of the boundary on the motion of a swimmer at 
a distance hz from the surface is 
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where C/ free and f2 Iree describe the swimmer’s free motion 
in an unbounded domain. For an axisymmetric spher¬ 
ical squirmer of radius a with a slip velocity u s (9) = 
J2n A nPn(0)r+J2n B nV n (9)G+J2 n C n V n (G)^>, as defined 
in [24] , this leading-order contribution to the motion due 
to this interaction is 
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where 9q is the angle the swimmer’s head makes with the 
wall normal and x is the direction of the head projected 
on the wall. All other components of the motion enter 
at higher order. Of these, particularly interesting is the 
rotation about the wall normal Q z , which leads to circular 
trajectories of particles swimming parallel to the wall. 
When 0Q = 7t/2 this rotation is given by 


in an easterly sense in the northern hemisphere but in 
a westerly sense in the southern hemisphere, so that a 
swimmer with this squirming mode qualitatively resem¬ 
bles bacteria like E. coli which have a counter-rotating 
head and tail. Indeed, such bacteria are known to swim 
in circles close to boundaries [4] . The direction of the cir¬ 
cling is set by the direction of the head-tail rotation and 
is right-handed if C 2 is positive. The radius of curvature 
is given by 160U llee h 4 /a 3 C 2 , indicating that this effect is 
strongly localised at boundaries, as reported experimen¬ 
tally [5]. 

The interaction of swimmers with planar surfaces has 
been studied using a variety of techniques, such as bound¬ 
ary integrals [25] or singularity approximations of the 
swimming action urn- The present approach using the 
reciprocal theorem reproduces all previous results, but 
it does so in a complementary fashion; the approxima¬ 
tion is in the integration kernel, and the description of 
the swimmer motion remains exact. Eqs. 0-0 are the 
leading order interactions for an arbitrary slip velocity, 
with no restriction on axisymmetry or shape. Successive 
improvements can be given by using more accurate forms 
of the stress tensor in the conjugate problem. 

As a second illustration, using exactly the same inte¬ 
gration kernel we can also find the motion of a passive 
particle propelled by the flow created by an active planar 
surface. For instance, this might correspond to transport 
by a carpet of cilia Esum, for which a simple appropriate 
slip velocity is the metachronal wave u s = Uq sin (kx — 
uit)x. Since the particle is no-slip, the reciprocal theorem 
reads U 4> + fl T = — J plane Uq sin (kx — u)t)x cr ■ z, so 
that its instantaneous velocity is given by 

U x = UQe~ kh {\ — kh) su\{kx — cot), 

U z = — Uoe~ kh khcos(kx — cot), (7) 

fly = ^Uoe~ kh ksin(kx — cot). 

The effect of the surface is exponentially confined to a 
layer of a thickness comparable to the wavelength of 
the metachronal wave.The motion of a swimmer close 
to such an active wall is given by the linear superposi¬ 
tion of eqs. 0.0 and (7j) , the boundary components of 
eq. 0. Depending on the wave’s direction, the deflection 
of a squirmer parallel to a wall [211 may be enhanced or 
suppressed, as illustrated in fig. Q. 

The velocity Q is precisely the flow field produced by 
a waving sheet, as found by G. I. Taylor in his seminal 
analysis [II . This makes sense because a small, spherical 
particle acts as a tracer which measures the local flow. 
Indeed, this simple example shows that the reciprocal 
theorem should be viewed as a solution for the entire 
swimmer hydrodynamics, and not just the motion. The 
fluid velocity at a point x is given by 
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where C 2 is the coefficient of the azimuthal slip velocity where <x is the stress tensor for the conjugate problem 
with angular dependence sin 29. Such a flow circulates where the particles Si are all force and torque free and a 
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FIG. 2: Flow field of the Stokes drag problem for two discs of 
unequal size, (a) Equal and opposite forces in the direction of 
the red arrows, (b) equal and opposite torques, clockwise on 
the upper disc. The black arrows show direction and relative 
magnitude of translational velocities and colour denotes flow 
speed in dimensionless units. 


point force 3> is applied at x. This greatly extends the 
scope of the reciprocal theorem from what was originally 
envisaged; not only does it give the swimmer motion, but 
a parallel analysis yields the full hydrodynamics. 

This use of the reciprocal theorem allows any swim¬ 
mer problem to be solved, provided the appropriate con¬ 
jugate stress tensor is known. Fortunately many exact 
and approximate solutions of Stokes flow problems for a 
variety of geometries exist in the literature, so that the 
relevant stress tensor is available in many cases. As we 
have seen, the general three-dimensional case is easily 
solved using an approximate integration kernel, giving 
good asymptotic results. However, this approximation is 
not appropriate in the near-field, where the separation 
becomes comparable to the swimmer size. We will now 
use a classic two-dimensional exact result for the stress 
tensor to shed light on these cases. 


III. INTEGRATION KERNEL IN TWO 
DIMENSIONS 

In two dimensions an exact solution for the Stokes drag 
on two discs is available, in the case where there are equal 
and opposite forces <f> and equal and opposite torques T 
acting on the discs [15] H25H50] . We state it in terms of 
fundamental flow singularities to connect it to the exist¬ 
ing literature on image systems. 

It is convenient to express the problem in the com¬ 
plex plane z = x + i y, admitting a formulation of fluid 
dynamics in terms of analytic functions m- In this do¬ 
main, we place two non-slip discs of radii ri, r 2 along the 
imaginary axis at z\ and Z 2 respectively. The flow about 
these discs can be expressed elegantly as being due to 
flow singularities at two points z = ±ii? inside the discs, 


where Sto, Rot, Str and Dip are, respectively, a stokeslet, 
rotlet, stresslet and dipole, each of unit strength and lo¬ 
cated at the position denoted by the subscript, defined 
as in [32]; and K \, K 2 set the disc radii and positions of 
their centres as 
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A more thorough discussion of this solution will appear 
elsewhere [351 . 

The case of two finite-sized discs is given by choosing 
ki,K 2 < 1. However, there are two other configurations 
for which eq. © holds: if re,; > 1 disc i is inverted and the 
fluid region is between two non-concentric nested circles, 
while if K'i = 1 then disc i becomes an infinite planar 
boundary [29]. In the latter case, shrinking the other 
disc to a point (k —¥ 0) we recover the well-known image 
systems for stokeslets and rotlets close to an infinite no- 
slip wall which, remarkably, are the same as those in three 
dimensions [23] . 

Although eq. © is exact, the restriction that there 
is no net force or torque means the reciprocal theorem 
can only be used to calculate the relative motions of the 
discs. This is easy to see from the singularity structure: 
the equal and opposite stokeslets and rotlets contribute 
forces and torques which cancel, while the second term in 
the rotlet opposes the couple induced by the point-force 
pair. This result explains why the existing solution for 
Stokes drag on a disc near a wall [33] does not exhibit 
the Stokes paradox: the wall exerts force and torque on 
the fluid equal and opposite to the dragging of the disc, 
regularising the asymptotic flow. The same observation 
is made in three dimensions by Blake & Chwang [23] . 

For completeness, we give the motions of the discs un¬ 
der this dragging, 
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These expressions demonstrate an exchange symmetry 
in the motion of the two discs: the force drives oppos¬ 
ing translation but co-rotation, while the torque gives 
opposing rotations and co-translation, with equal mag¬ 
nitude when the discs are of equal size. They also show 
explicitly that as the disc radii diverge they become im¬ 
mobile. 
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IV. APPLICATION OF THE SOLUTION 

Given the exact solution to the Stokes drag, eq. (|9), it 
is straightforward to construct an exact integration ker¬ 
nel to determine the interaction of any pair of active or 
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FIG. 3: Flow and trajectories of a finite-sized squirmer due to 
an active wall in two dimensions. In the absence of wall activ¬ 
ity the squirmer is deflected (black line). If the wall generates 
a metachronal wave, here with a wavelength comparable to 
the squirmer’s size, the deflection is enhanced or suppressed 
respectively when the wave opposes (blue) or co-moves (red) 
with the swimmer. This is qualitatively similar to the three- 
dimensional approximate case shown in fig. [I] 


inactive circular boundaries, for any slip velocities. The 
only limitation is that only the relative motion may be de¬ 
termined, since we cannot specify the forces and torques 
independently. If, however, one of the two boundaries 
corresponds to an immobile object the behaviour of the 
other may be determined completely. 

Consider first a passive disc of radius a advected by an 
active wall (k\ —> 1) with a surface metachronal wave, as 
before. The motion of the disc is given by U ■ + f IT = 
— f_ dx Uq sin (kx — ut)x ■ cr ■ z, or, explicitly, 


Ret/ = U 0 e~ kVh2 - a ' 2 sin (kx - tut), 

Im U = — Uoe~ k ^ h2 ~ a2 k(h — a 2 /h) cos (kx — ut), (14) 
ft = U 0 e~ k ^ h2 ~ a2 1 — a 2 /h 2 sin (kx — ut), 

and is broadly the same as eq. ([TJ) in the three- 
dimensional case. The finite size of the tracer is encoded 
by the replacement of h with \Jb? — a 2 . When a-lOwe 
recover, as before, an exact flow field external to the wall 
due to the activity, shown in fig. [3| 

For a squirming particle we adopt a minimal model 
consisting of the two squirming modes determining the 
swimming speed V and the force dipole W, so that the 
slip velocity is given by 

u s = 2V sin[0 - 0] + 2W sin [2 (<j> - 0)], (15) 


at the point on the surface parametrised by the angle <f> 
with respect to the direction of the head, 0. When V 
and W have the same sign the swimmer is termed con¬ 
tractile, otherwise, it is extensile [14] : these two types of 
activity are related by time inversion [34] . Applying the 
reciprocal theorem for this slip velocity gives the motion 



FIG. 4: A squirmer of radius V 2 = 0.1 in a tank of radius 
ri = l. (a) An extensile swimmer rapidly adheres to the wall, 
(b) A contractile swimmer with the same initial orientation 
and squirming strengths takes a longer time but also settles on 
a trajectory following the boundary in the opposite sense to 
the extensile swimmer, (c) The gap between the swimmer and 
the wall for the two trajectories. The asymptotic separation 
is significantly larger for the contractile swimmer. 


as 
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We give this expression for the general case of two ar¬ 
bitrary circular boundaries with slip velocities given by 
eq. ©, from which specific cases can be recovered by 
appropriate choices of K \, k 2 - Letting aci —>■ 1, and 
V], W\ 0, gives the interaction of a squirming disc with 
a wall [18]. Superposition with eq. (14) gives the motion 


of a squirmer in the presence of a metachronal surface 
wave. The general phenomenology in this case is illus¬ 
trated in fig. [ 3 ] when swimming against the metachronal 
wave the deflection from the surface is increased (blue 
line), while swimming with the metachronal wave de¬ 
creases the deflection (red line). Qualitatively, this is 
the same as we found in three dimensions. 

Next, consider a swimmer inside a no-slip circular tank. 
If V = 0 the swimmer is apolar and we recover previous 
results [35] (not shown). In particular, the symmetry 
of the problem means such a swimmer gets stuck in the 
centre of the tank. Including self-propulsion, the swim¬ 
mer crashes into the boundary and settles into a trajec¬ 
tory following it, with the reorientation and sense of the 
trajectory determined by the microscopic details of the 
squirming, as shown in fig. [4] For the initial conditions 
studied, extensile swimmers are found to be attracted 
more strongly to the boundary than contractile ones. 
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FIG. 5: (a) Phase flow for a swimmer with radius 0.5 relative 
to the post, (b) Phase flow for a swimmer with equal radius 
to the post. When the swimmer radius is smaller than the 
post radius there are closed trajectories about a circular orbit 
(blue). As the swimmer radius approaches the post radius 
from below the distance of this circular orbit diverges as (1 — 
P 2 /Vi) _1 ^ 2 ■ The exclusion is shown in green, (c) Existence 
of the circular orbit (blue) compared to the excluded radius 
(green). 


We now turn to interactions between two finite-sized 
squirmers. The fact that we can only find relative trans¬ 
lation and rotation becomes problematic in determining 
the motion; in the previous cases a frame of reference was 
defined by the wall or tank. There are two resolutions: 
to make one disc an immobile, passive ‘post’ which fixes 
a reference frame as before; or to make both squirmers 
identical so that each contributes equally to the motion in 
some centre-of-mass frame, which cannot be determined. 


In the former case, taking Vj, W\ —»• 0 disc 1 becomes 
a no-slip post. Then we can transform eqs. (16)-(18) to 
polar coordinates (d, ip) centred on the post. By the rota¬ 
tional symmetry about the post, this dynamical system 
depends only on the relative angle Or = 9 — <p. This al¬ 
lows us to form a lower-dimensional system, (d(f), 0n(t )), 
integrable by separation of variables. A fixed point in 
this phase space giving a bound trajectory has previously 
been seen for the interaction of a swimmer with an infi¬ 
nite wall m G23 ■ We find that this effect persists when 
the post is finite in size, corresponding to a bound orbit 
of the post by the swimmer. As the size of the swimmer 
approaches the size of the post, the location of this orbit 
diverges algebraically as d ~ (1 — r^/ri) -1 ' 2 ■ When the 
swimmer is larger than the post we see generic deflection. 


Finally, consider the hydrodynamic interaction be¬ 
tween two identical squirmers and an initial condition 
with a mirror symmetry. Then each will rotate in an op¬ 
posite sense with respect to the centre of mass, with half 
the magnitude of fl 2 ~ fii- We illustrate in fig. [6] some of 
the close proximity interactions, determined by exact hy¬ 
drodynamics. The reorientation and distance of closest 
approach for a typical collision is shown in fig.[6j(a) , where 
the squirmers start exactly parallel. Angling two con¬ 
tractile swimmers slightly towards each other can lead to 
bound trajectories, shown in blue in fig.[6](b). This occurs 
for extensile swimmers when they initially point away 
from each other, as a consequence of the time-reversal 
symmetry between the two types ]3?J- 



FIG. 6: Interactions of two identical squirmers in the centre- 
of-mass frame, swept out in time, (a) Two initially parallel 
extensile swimmers are attracted and undergo the reorienta¬ 
tion shown by the red trace. The minimum separation during 
the collision is shown by the grey curve, (b) Trajectory of 
two extensile (red) and contractile (blue) swimmers with an 
initial tilt towards each other. Contractile swimmers exhibit 
a periodic bound trajectory. By the time-reversal symmetry 
of contractile and extensile activity, this is seen in extensile 
swimmers when the initial orientation is outwards. 


V. DISCUSSION 

The reciprocal theorem is a standard technique for un¬ 
derstanding swimmer motion |l5W20j . The reinterpreta¬ 
tion presented here shows how it can be also used to 
obtain all aspects of swimmer hydrodynamics. We gave 
simple illustrations in three dimensions for the motion of 
a squirmer in the presence of a planar boundary, repro¬ 
ducing literature results without approximation of the 
swimming strategy, and extending to interactions with 
a ciliated surface. Adopting standard exact results for 
Stokes flows in two dimensions [3S] allowed for a wide 
range of swimmer hydrodynamics and near-field inter¬ 
actions to be solved. The three-dimensional results are 
in good qualitative agreement with these exact results, 
suggesting that one may have some confidence in their 
validity, even close to surfaces. 

If the stress tensor of a Stokes drag problem is known, 
the corresponding swimmer hydrodynamics is also given 
by the reciprocal theorem. This allows a large number of 
swimmer hydrodynamics problems to be solved exactly in 
both two and three dimensions, by adoption of literature 
results [33]. Another interesting generalisation would be 
to incorporate free surfaces with boundary conditions of 
specified stress rather than slip velocity mm- 

The reciprocal theorem gives swimmer hydrodynamics 
in terms of an integral against the boundary data, with 
the integration kernel identified as the stress tensor of a 
conjugate Stokes drag problem. This is analogous to, but 
distinct from, boundary integral methods [38] for hydro¬ 
dynamics, where the flow is given by integrating a force 
distribution over the boundaries. Here, the boundary 
conditions are the slip velocities directly, and the integra¬ 
tion kernel is the stress tensor for a conjugate Stokes drag 
problem. Because the stress tensor decays more rapidly 
than the Green function for a point force, this approach 
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